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I .  INTRODUCTION 


The  effect  of  particle  motion  on  spectral  lineshape  (i.e.t  Doppler 
broadening)  is  usually  evaluated  by  assuming  that  the  particles  have  a  random 
thermal  motion.^  In  high-speed  gas-flow  lasers,  of  which  chemical  lasers  are 
an  example,  the  working  fluid  may  have  significant  mean  velocities  as  well  as 
random  velocities  in  the  optical  path  direction.  In  these  cases,  it  is  neces¬ 
sary  to  take  the  mean  motion  into  account  when  evaluating  the  spectral  line- 
shape.  The  effect  of  a  source  flow  on  Doppler-broadened  lineshape  is  eval¬ 
uated  herein  for  the  case  where  radiation  is  perpendicular  to  the  source  flow 


axis 


I  — 


The  application  to  a  cw  chemical  laser  is  then  noted. 


f 


*A.  Yariv,  Introduction  to  Optical  Electronics,  (Holt,  Rinehart  and  Winston, 
New  York,  19^1)  pp.  76-79. 

^D.  J.  Spencer,  H.  Mirels,  and  D.  A.  Durran,  J.  Appl.  Phys.  43,  1151  (1972). 
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II.  THEORY 


Consider  a  two-dimensional  source  flow  with  half  angle  0g  (Fig.  1).  A 
radiation  field  I(v)  is  assumed  to  propagate  in  a  direction  that  is  transverse 
to  the  flow  axis.  The  divergence  of  the  flow  results  in  mean  motion  in  the 
transverse  (optical  path)  direction.  Local  mean  velocity  and  random  velocity 
in  the  transverse  direction  are  denoted  by  V  and  v,  respectively.  The  resul¬ 
tant  transverse  velocity  field  is  W  =  V  +  v.  The  radiation  frequency  v,  which 
is  resonant  with  particles  of  velocity  W,  is  found  from  the  Doppler  relation^ 

W  -  (v  -  vq)  c/vo  (1) 

where  vq  is  the  resonant  frequency  for  stationary  particles  and  c  is  the  speed 
of  light.  The  spectral  lineshape  is  found  by  determining  the  distribution 
function  for  particles  in  the  range  W  to  W  +  dW,  which  is  found  in  the  fol¬ 
lowing  paragraphs. 

2 

Consider  the  case  0g  «  1.  The  flow  density  and  axial  flow  velocity  at 
each  streamwise  station  x  are  independent  of  0,  whereas  the  transverse 
velocity  can  be  expressed 


(2) 


where  Ve  is  the  value  of  V  corresponding  to  9  *  0e-  The  fraction  of  particles 
in  the  velocity  range  V  to  V  +  dV  is 


F(V)  -  1/(2  Ve)  IVl  <  Ve 

-  o  I vl  >  ve  {3) 

Under  equilibrium  conditions,  the  fraction  of  particles  in  the  velocity  range 
v  to  v  +  dv  has  a  Maxwellian  distribution* 

i  2.2 

f(v>  -  -i-  «-V  /a  (4) 

ir  a 

where  a  =  (2  kT/m)*^  is  the  most  probable  random  particle  speed.  The 
quantities  F(V)  and  f(v)  are  normalized  so  that 


£ 

£ 


f(v)dv  =  r°  F(V)  dv  =  1 


f ( v)  F(V)dv  dV 


1 


(5a) 


(5b) 


The  integrand  in  Eq.  (5b)  denotes  the  fraction  of  particles  in  the  combined 
range  V  to  V  +  dV  and  v  to  v  +  dv.  We  now  replace  the  independent  variables 
V,  v  by  the  pair  W,  w  defined  by 


W  -  v  +  V 


(6a) 


8 


w 


v  -  V 


(bb) 


It  is  seen  from  Fig.  2  that  W,  w  are  orthogonal  variables.  The  particles  are 

confined  to  the  region  -V  <  V  <  V  in  V,  v  space  and  to  the  region  W  -  2  V 

<  w  <  W  +  2  Vg  in  W,  w  space.  The  particles  in  the  velocity  range  V  to  V  +  dV 

and  v  to  v  +  dv  are  related  to  the  particles  in  the  range  W  to  W  +  dW  and  w  to 

w  +  dw  by 


f(v)  F(V)  dvdV  =  f(— -  )  F(~) 


a(y,v) 

3(w,W) 


dw  dW 


(7) 


o 

where  3(v, V)/3(w,W)  is  the  Jacobian  given  by 


3(v,V)  _  iVw  VWi  _  1 
3(w,W)  “  vw  Vw‘  =  2 


Equation  (5b)  becomes 


/*  /*W+2 

dW  J 

—  W-2V 


W+2Ve  -(vrW)2/(4a2) 

dw  (£ - m - 

4ir  '  aV 


) 


£ du  [ 


w+v  w-v 

erf( - e)  -  erf( - e) 

v  a  1  v  a  '  i 

4V 


(8) 


(9a) 


(9b) 


The  integrand  in  Eq.  (9b)  can  be  interpreted  as  the  fraction  of  particles  in 
the  range  W  to  W  +  dW.  Introduce  the  notation  Xg  =  Vg/a,  and 

•^D.  H.  Menzel,  Fundamental  Formulas  of  Physics,  Vol.  I,  (Dover  Publications, 
New  York,  1960) ,  p.  26^ 
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(10) 


,,  v-v  ,  /0  v-v 

X  =  H  =  - ®  -c-  =  2 (in  2)1/2  ° 

a  v  a  v  ' 


Av„ 


where  Av^  Is  the  Doppler  width  (FWHM)  corresponding  to  Xg  =  0.  The  fraction 
of  particles  in  the  interval  X  to  X  +  dX  is  denoted  n(X,  Xg)  and,  from  Eq. 
(9b),  equals 


n(X,  X  )  =  [erf  (X  +  Xe)  -  erf(X  -  Xe)]/(4X  )  (11) 


Limiting  forms  of  n(X,Xg)  are 


(12a) 


-  e'X  [1  +  0(X2)] 


(12b) 


,l/2 

2X 


erf  (Xe)[l  +  0(X2)j 


(12c) 


2 

Thus,  transverse  mean  flow  effects  are  negligible  for  (Xg)  <<  1  and  must  be 
considered  for  (Xfi)2  >  0(1).  This  result  is  physically  realistic,  since  Xe  is 
the  ratio  of  characteristic  mean  to  characteristic  random  motion  in  the  trans¬ 
verse  direction.  The  variation  of  *1/2  n(X,Xe)  with  X  is  plotted  in  Fig.  3  for 
various  values  of  Xg.  Figure  3  can  be  used  directly  to  estimate  the  effect  of 


11 


transverse  mean  motion  on  zero  power  gain  of  a  Doppler-broadened  medium.  If 


we  let  gQ(X,  Xe)  denote  the  zero  power  gain  corresponding  to  X,  Xe,  It  follows 

that  In  the  Doppler  limit4*5  (l.e..  In  the  limit  of  Av /Av  <  <  1  where  Av. 

n  D  n 

Is  the  homogeneous  half-width) 


g0(x,xe) 

g'o'tO'.f)) 


X1'2 


n(x,Xe) 


(13) 


where  go(0,0)  Is  the  line  center  value  In  the  absence  of  transverse  motion. 
Thus,  the  ordinate  In  Fig.  3  Is  a  direct  measure  of  the  effect  of  transverse 
mean  motion  on  zero  power  gain.  Line  center  gain  is  decreased,  and  the 
lineshape  Is  broadened  as  X_  is  increased.  The  corresponding  value  of  the 
anomalous  index  of  refraction  n(X,Xfi)  can  be  expressed,4,5  for  X  =  C/VQ» 


2. 

X  ■  g(0,0) 


dX 

rf  [erf(Xo+Xe)  -erf(Xo-xJ] 

O 


■~i72  r  D(Vt‘  +  0(xZ)! 


(14a) 


(14b) 


“T72  D(X)[1  +  0(X^)]  (14c) 


V  Mi re Is,  AIAA  J.  (5),  478  (1979) 

5M.  Sargent,  M.  0.  Scully  II,  and  W.  E.  Lamb,  Jr.,  Laser  Physics,  (Addlson- 
Wesley,  Massachusetts,  1974),  pp.  144-155. 
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where  D(  )  is  the  Dawson  integral.  Equation  (14a)  is  plotted  in  Fig.  4.  The 
maximum  value  of  the  index  decreases  as  Xg  increases  and  occurs  at  values  of  X 
and  Xe  related  by  D(X  +  Xe)  -  D(X  -  Xe),  which  for  (X  +  Xe)“2  «  1  becomes  X2 

=  y}  +  0.5  [i  +  o(x  +  xer2]. 
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III.  APPLICATION 


It  is  convenient  to  express  the  parameter  Xg  in  terms  of  local  flow  Mach 
number.  The  most  probable  speed  a  is  related  to  the  local  speed  of  sound  a8 
-  (yRT/m)1^2  by  the  relation  ag/a  ■  (y/2)*^2  where  y  is  the  ratio  of  specific 
heats,  which  equals  7/5  and  5/3  for  diatomic  and  monatomic  gases,  respectively. 
Let  Mg  5  Ve/ag  denote  the  characteristic  transverse  flow  Mach  number,  and  it 
follows  that  Xe  is  related  to  Me  by 


xe  -  (2/y)1/2  Me  (15) 

Thus,  Xe  is  nearly  equal  to  the  characteristic  transverse  flow  Mach  number. 

The  present  results  can  be  applied  to  cw  chemical  lasers  in  which  trans- 
verse  flow  expansion  is  permitted.  Typical  axial  flow  Mach  numbers^  are  in 
the  range  4  to  6  and  can  be  characterized  as  being  of  order  5.  The  cor¬ 
responding  value  of  Xe  is  then  of  order  Xe  =  0(50^).  Thus,  transverse  flow 

effects  are  negligible  for  02  «  0(1/25)  and  must  be  considered  for  0^  >  0(1/25). 

e  e 

The  effect  of  source  flow  on  cw  chemical  laser  output  power  is  being  inves¬ 
tigated  for  the  case  of  a  two-level  model, ^  laminar  mixing,  a  Fabry  Perot 
resonator,  and  multiple  longitudinal  modes  spacing  [i.e.,  Av^/Av^  <0  (1), 
where  Av^  ■  c/2L  is  the  longitudinal  mode  spacing].  Preliminary  results  are 
indicated  in  Fig.  5.  The  ordinate  P/Psat  is  the  ratio  of  output  power  to 
saturated  output  power,  and  the  abscissa  is  the  ratio  of  cavity  threshold 
gain,  gc,  to  the  maximum  zero  power  gain,  g^  ,  for  the  case  Xg  «  0.  The 
decrement  in  output  power  caused  by  source  flow  is  seen  to  become  more  severe 
as  threshold  gain  is  increased. 
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